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Abstract 

We consider a Ginzburg-Landau type energy with a piecewise constant pinning term a in 
the potential (a^ — |itp)^. The function a is different from 1 only on finitely many disjoint 
domains, called the pinning domains. These pinning domains model small impurities in a 
homogeneous superconductor and shrink to single points in the limit e — > 0; here, e is the 
inverse of the Ginzburg-Landau parameter. We study the energy minimization in a smooth 
simply connected domain Q C C with Dirichlet boundary condition g on dfl, with topological 
degree degg^^g) = > 0. Our main result is that, for small e, minimizers have d distinct 
zeros (vortices) which are inside the pinning domains and they have a degree equal to 1. The 
question of finding the locations of the pinning domains with vortices is reduced to a discrete 
minimization problem for a finite-dimensional functional of renormalized energy. We also find 
the position of the vortices inside the pinning domains and show that, asymptotically, this 
position is determined by local renormalized energy which does not depend on the external 
boundary conditions. 
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1 Introduction and main results 

In this work we study the minimizers of the Ginzburg-Landau type functional 

EsA^) = mi^^i' + ^(«' - 1^1')'} ' 

where C C is a bounded, smooth, simply connected domain, e is a positive parameter (the in- 
verse of the Ginzburg-Landau parameter k = 1/e), S = 6{e) > is a geometric parameter and u 
is a complex- valued map. In order to define the function ag, we need to introduce the notion of a 
pinning domain. 

Fix M € N* points ai, um £ ^- Let uj be an open subset such that uJ C B(0, 1) and G w. 
For 1 < i < M and for all 5 > denote := + S ■ uj, i.e. the set to scaled by 6 and centered at 

Definition. The set cos := ufliCOg is called a pinning domain. 

- — I \M nl,, _ 



For example, \i uj = B(0, ^), then the pinning domain is cos = U^£]^i?(aj 



We now define ag : 0, ^ {b,l}, b £ (0, 1) as: 

II it X £ \l \uJs 

The functionals of this type arise in models of superconductivity for composite superconductors. 
The experimental pictures suggest nearly 2D structure of parallel vortex tubes (|25). Fig 1.4). There- 
fore, the domain can be viewed as a cross-section of a multifilamentary wire with a number of 
thin superconducting filaments. Such multifilamentary wires are widely used in industry, including 
magnetic energy-storing devices, transformers and power generators |17) . |16) . 

Another important practical issue in modeling superconductivity is to decrease the energy dissi- 
pation in superconductors. Here, the dissipation occurs due to currents associated with the motion 
of vortices (t21j, |i6j). This dissipation as well the thermomagnetic stability can be improved by 
pinning ("fixing the positions") of vortices. This, in turn, can be done by introducing impurities or 
inclusions in the superconductor. In the functional ([T|) the set ojs models the set of small impurities 
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in a homogeneous superconductor. The size of the impurities in our model is characterized by the 
geometric parameter 5 which goes to zero together with the material parameter e. We assume 
henceforth that 



Essentially, this condition means that 5 is much larger than e on the logarithmic scale. For example, 

if e = 2-J and 6(e) = 2-^(-'), then (|H]) implies that 0. 

j 

Notation. In what follows: 

• We consider a sequence i and we write e instead of the dependence of e on n is implicit. 

• We simply write 5 (instead of (5(e)); the dependence of 5 on e is implicit. 
We study the minimization problem for the functional ([1]) in the class 



where g £ C°°(0r2,S^) is such that deggQ{g) = d > 0. Recall that the degree (winding number) of 
g is defined as 



Here "x" stands for the vectorial product in C, i.e. zi x Z2 = Im{'ziZ2), zi,Z2 G C, and dr is the 
tangential derivative. The degree is an integer, and the condition degQ^{u) = d > 0, u G H^{Q,C) 
implies that u must have at least d zeros (counting multiplicity) inside Q. The properties of the 
topological degree can be found, e.g., in |13) or [8]. 

Minimization problems for Ginzburg-Landau type functionals have been extensively studied 
by a variety of authors. The pioneering work on modeling Ginzburg-Landau vortices is the work 
of Bethuel, Brezis and Helein In this work the authors suggested to consider a simplified 

Ginzburg-Landau model ([T|) with a = 1 in Q (ie. without pinning term), in which the physical 
source of vortices, the external magnetic field, is modeled via a Dirichlet boundary condition with a 
positive degree on the boundary (|3]). The analysis of full Ginzburg-Landau functional, with induced 
and applied magnetic fields, was later performed by Sandier and Serfaty in |27) . 

The functional ([1]) with non-constant a{x) was proposed by Rubinstein in |26j as a model of 
pinning vortices for Ginzburg-Landau minimizers. Shortly after, Andre and Shafrir [3] studied the 
asymptotics of minimizers for a smooth (say C^) a. One of the first works to consider a discontinuous 
pinning term, which models a composite two-phase superconductor, was [18) . In this work, a single 
inclusion described by a pinning term independent of the parameter e was considered for a simplified 
Ginzburg-Landau functional with Dirichlet boundary condition g on d^l. Namely the pinning term 
is 



here w is a simply connected open set s.t. tJ C fi. The main objective of |18) was to establish that 
the vortices are attracted (pinned) by the inclusion uj, and their location inside oj can be obtained 
via minimization of certain finite-dimensional functional of renormalized energy. Full Ginzburg- 
Landau model with discontinuous pinning term was later considered by Aydi and Kachmar [5]. An 
e-dependent but continuous pinning term as{x) was studied by Aftalion, Sandier and Serfaty in [T]. 
The work [3j studies the case of the smooth a with finite number of isolated zeros, and in [2j the 
pinning term a takes negative values in some regions of the domain fi. The other works related to 
Ginzburg-Landau functional with pinning term include, e.g., |21) . |28) . 




(H) 



Hi := {u€H\n,C)\tT9nU = g} 



(3) 
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In this work, we consider the minimization problem ([I])-® with a discontinuous pinning term 
given by (|2]). We prove that despite the fact that — )• 1 a.e. as e — )• 0, i.e. the pinning term 
disappears in the hmit, the pinning domains us capture the vortices of Ginzburg-Landau minimizers 
of ([1]) for small e. 

The main difficulty in the analysis of this problem stems in the fact that the a priori Pohozhaev 
type estimate ||1 — |'1'P||^2(q) < Ce"^ for the minimizer v (on which the analysis in and [18] 
is based) does not hold. Therefore, we develop a different strategy of reducing the study of the 
minimizers of ([!]) to the analysis of §^-valued maps via the uniform estimates on the modulus of 
minimizers away from the pinning domains (see Proposition |5] below) . 

Following |18j . let Us be the unique global minimizer of in with f/g = 1 on dO,. This C/g 
satisfies b < < 1. For v G Hg we define 

Using the Substitution Lemma of |18j . we have that for v G Hg^ 

E,{UeV)=E,{Ue)+F,{v). (4) 

From the decomposition (|4|), we can reduce the minimization problem ([I])-® to the minimization 
problem for Ffr in H^., namely, the minimizer of in Hg has the same vorticity structure as the 
original minimizer of ([H)-©. 

Depending on the relation between M (number of inclusions), and d (number of vortices), we 
distinguish two cases: 

Case I: M > d (more inclusions than vortices). 
Case II: M < d (more vortices than inclusions). 

For example, we are going to show that for the minimizer Ve'- 

• if M = 3 and d = 2 (Case I), we have two distinct inclusions containing exactly one zero each, 

• if M = 2 and d = 2> (Case II), we have one zero inside one inclusion and two distinct zeros inside 
the other inclusion. 

Generally speaking, outside a fixed neighborhood of d' = min {d, M} inclusions (centered at a = 
(oj^, Oj^, )), the minimizer Ve is almost an S^-valued map. Moreover, by minimality of ti^, the 
selection of centers of inclusion containing its zeros and the distribution of degrees of around the 
Cj's are related to the minimization of the Bethuel-Brezis-Helein renormalized energy Wg. In other 
words, we reduce the problem of finding vortices of the minimizers Vf, to a two-step procedure. As 
the first step, we determine the inclusions with vortices, which is a discrete minimization problem 
for Wg and is significantly simpler then the minimization of this renormalized energy functional 
over 0"^'. Secondly, we determine the locations of the zeros (vortices) locally inside each inclusion 
and show that their positions depend only on 6, on the geometry of oj and on the relation between 
d and M, but not on the external Dirichlet boundary condition g (see Theorem H] below) . 
Our main result in Case I is the following: 

Theorem 1. Assume that M > d. Let he a minimizer of F^ in Hg{^). For any sequence Sn i 0, 
possibly after passing to a subsequence, there are d distinct points {oj^, ...,aj^} C {cj, 1 < i < M} 
and a function v* G H^^^{^ \ {ojj, ...,ai^},§^) such that: 

1. V* is a harmonic map, i.e. 

{-Av* = v*\\7v*\'^ m Q\{ai-^, ...,ai^} 
V* = g on do, 
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(5) 



2. Wehavevs„^v* strongly in Hl^^{Q.\{ai^, ...,ai^]) and Ve„ ^ v* m CgJ.(ri \ {ai, om})- 

3. Vs„ has d distinct vortices such that is inside w^™, m = l,...,d and for small 
fixed p, degQB(^^n p^{ve„) = 1. 

^. The following expansion holds 

F,{ve) = ndb^l lne| + tt{1 - b^)d\ ln5\ + Wg{{a^,^), (a^,, 1)) + W + Oeil). (6) 

Here W > is a local renormalized energy depending only on d, b and cj . Moreover, the 
d-subset {cjj, ...,0^^} C {ai, ...jOm} minimizes the Bethuel-Brezis-Helein renormalized energy 
Wg among the d-suhsets o/ {oi, om}- 

Remark 1. Here, Wg denotes the renormalized energy given by Theorem 1.7 in [IjJ (with the degrees 
equal to 1 and the boundary data g). Its definition is recalled in Section [5.31 

The main result in Case II is 

Theorem 2. Assume that M < d. Let be a minimizer of in Hg{Q). For any sequence En i 0, 
possibly after passing to a subsequence, there is v* G //^^^(ri \ {ai, oa/}, S^) which satisfies ([5]) in 
O \ {ai, Om}; such that: 

1. Ve^^v* strongly m Hl^^{Vt\{ai,...,aM}) and v^^ ^ v* in Cj^ (il \ {oi, aA/})- 

2. For p > small, v^^ has exactly di := <i&SdB{ai,p)i'^e„) zeros in B{ai,p). They are isolated, lie 
inside cj^ and they have a degree equal to 1. 

3. 

, \ d ] 7 r ^ 1 7 r d 

< di < + 1, where — is the integer part of —. (7) 



M 



M 



M 



Moreover, if -§j = 'mo G N, then di = mo, 1 < i < M . Otherwise, the configuration 
{(ai, di), ...,(aAf , dAf)} minimizes the renormalized energy Wg among the configurations 
{(ai, di), (oAf , dA/)}- Here {aj \1 < i < M} are fixed and di € Z are the subjects to the 
constraints ^ and ^fl^ di = d. 

4- The following expansion holds when e — )■ 

M 

ini Fs = TTdb^\lne\ + TT{^d^, - db'^)\ln6\ + Wg {{a, d}) + W + Oe{l). (8) 

^9 i=l 

Here, {a,d} = {(ai, di), (oAf, c^Af )} is a configuration given by the previous assertion and 
W is local renormalized energy which depends only on uj,b,d and M. 

In both cases, we prove that the asymptotic location of the vortices inside a pinning domain 
depends only on 6, cj and on the number of zeros inside the inclusion (see Theorem |4]) : this location 
is independent of the boundary data g on OO. 

2 Main tools 

In this section we establish: 

• Estimates for f/^, 

• Upper bounds for the energy of minimizers in Case I and Case II, 

• An ry-ellipticity estimate for minimizers. 
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2.1 Properties of 

Proposition 1 (Maximum principle for Ue, [18\ Proposition 1). The special solution Us satisfies 
b<Ue <l inn. 

Proposition 2. There are C, c > (independent of e) s.t. for any R > we have 

\as - Ue\ < Ce""^ in Vr:={x eQ] dist(x, dujs) > R}, (9) 

\VUe\ < — - — m Vr. (10) 
The proof of the Proposition [5] is presented in the Appendix |X] 

2.2 Upper Bounds 

Proposition 3. Let ^ = ^■ 



1. Upper hound in Case /; M > d 

There is a constant C depending only on g,uj and Q s.t. we have 

inf FJ-,n) <7rdb^\lnC\+TTd\ln6\+C. (11) 

2. Upper hound in Case II: M < d 

There is a constant C depending only on g,u] and i7 s.t. for all di, ...,dM S N s.t. Y^di = d 
we have 

inf Fe{-,n) <Trdb'^\lnC\+Tr^dj\ln6\+C. (12) 
The proof of Proposition [3] is given in Appendix [B] 

2.3 Identifying bad discs 

Lemma 1. Let gs,go G C°°{dn,C) he s.t. < 1 — l^f^l < e and g^ — )• go in (7^(50). Let also 
as,l3s€L°-{n,[b, 1]). 

Consider the weighted Ginzhurg-Landau functional 

i^rw = ^/J«.iv.p + |(i-H2)2}. 

Denote a minimizer of in H^^ . Then the following results hold: 

1. Let X = Xe G (Oi 1) be s.t. x ~^ 0. There are eo > 0, C > and Ci > depending only on 
b,X,^,\\9o\\cHdn) s.tfore<eo, if 

F^ive,B{x,e'/^) nn) < x^\lne\ - Ci, 

then 

\vs\ >l-Cx m B{x,e^/^)nn. 

2. Let fi G (0,1). Then there are eo,C > depending only on b,iJ,,n, \\go\\c'^{dn) ^-t- for e < Sq, 

F^{vs,B{x,e^^^)nn) < C\lne\, 

then 

\vs\ > fi in B{x,£^/'^) n 9.. 

Lemma [1] is proved in Appendix O 



3 A model problem: one inclusion 

By combining the results of Section |2l the proofs of both Theorem [1] and Theorem [2] are based 
on the analysis of two distinct problems: 

1. A minimization problem of the Dirichlet functional among §^-valued map defined on a perfo- 
rated domain. 

2. The study of the minimizers around an inclusion. 

This section focuses on the second problem. More precisely, we fix /) > and study the minimization 
problem of F£{-, B{ai, p)) with variable boundary conditions. 

Fix p > and let fe, fo £ C°°{dB{0,p)) be s.t. /o is S^-valued and s.t. 

Wfe - /o||ci(9B(0,p)) (13) 

and 

\\\fe\-l\\L^iaBio,p))<Ce\ (14) 

Assume also that deggs^o^p^ife) = degg£(o,p)(/o) = do > 0. 
For i G {1, ...,M} consider the minimization problem 

F,{v,B{a„p)):=^ j iu^^\Vv\^ + ±^U^{l-\v\^f\Ax (15) 

in the class 

Hl^, := {v G H\B{ai,p)X)\tTQBia,,p)v{x) = /.{x - a^)}. (16) 

Without loss of generality assume Oj = 0. Let be a minimizer of (|15p in (|16p . Performing the 
change of variables x = | in (|15p . we have 

F,{ve, B{0, p)) = F^{v,, BiO, ^)) := ^ J^^^ S^U!\Vv\' + ^U^l " l^l')'} dx. (17) 
Here, for a map w S H^{B{0, p)), we denote w{x) := w{5x) and ^ = -. The class (|16p under this 



change of variables becomes 

Hi := [v G H\B{^, ^), C) | trgB(o,£)7)(-) = /.(5-)} . (18) 

Note that the above rescaling enables us to fix the pinning domain independently of e. 
The asymptotic behavior of will be obtained in several steps: 

• We first establish a bound for \Vf,\ (Proposition E)). This bound will allow us to localize (roughly) 
the vortices of near the inclusion. 

• We next establish sharp energy estimates (Proposition [6| and use them to obtain the uniform 
convergence of solutions away from the inclusion (Proposition [7] and Corollary |2]) . We establish 
the strong convergence of solutions away from the "vortices" (Proposition [8| and derive the 
equation satisfied by the limiting map (Proposition [TO]) . 

• The last step is the location and quantization of the vorticity: for small e, the minimizers admits 
exactly cLq zeros, and all the zeros lie in the inclusion and have a degree equal to 1 (Propositions 
[8]and[II]). 

Following the same lines as for Proposition [3l one may prove 

Proposition 4. Let Vf, he a minimizer of F^ in il6\) . Then there is a constant C independent of e 
s.t. we have 

Fe{ve,B{0,p)) = Fi:{v„B{0, ^)) < ndob^l In^l + ird^ ln6\ + C. (19) 
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3.1 Uniform convergence of Iv^l to 1 away from inclusions 

Proposition 5. Let K CM."^ be a compact set such that u d K and dist(5X, cj) > 0. Then there is 
C > independent of e s.t. for sufficiently small e we have 

\ve\ > 1-C|lner^/^ m Be\K. 

Proof. Using Lemma [1] with x = I Inej""'^^'^, we find that there exist C, Ci > s.t. for e > small, 
if F,{ve,B{x,e^/^)) < |lne|^ - Ci then \ve\ > 1 - Cx in B{x,e^/^). 

We argue by contradiction. Assume that there exists a compact K containing uj s.t. dist((9i^', lo) > 
and s.t., up to a subsequence, there is a sequence of points G B{0, ^) \ K s.t. < 1 — 

C\ lne|~-^/^ with C given by Lemma[T] Note that Xe £ B{0, ^)\K corresponds to Xe G B{0, p)\{5-K). 
From Lemma [1] and Proposition |2] 



1 [ ||V.,|2 + -i^(l-kn4>|lne|V3_0(l). 

2 yB(x„ei/4) t 2e^ J 



(20) 



We claim that due to the conditions (|13p . (|14p . we may extend (keeping the same notation for 
the extension) to a smooth map, still denoted v,., s.t. 

've{x) = x'^o/\x\do in B{0,3p) \ B{0,2p) 

(1 - \v,\^f < Ce^ . (21) 



S(0,3p)\B{0,p) 

Vi^el < C with C > is independent of e 

To make the above extension explicit, choose Q G C°°(M^, [0,1]) s.t. = in [0,p] and Q = 1 m 
[2/9, 3/9] and take 



v,{sn= C(.) + (l-C(s))l/.(pe^')l 



^*[doe+(l-C(s))<^.(pc»»)] 



Here x = se'^ , s > and (/.^ G C°°(aB(0, />), M) is s.t. fe{pe'^) = \ f^\e'('ioe+4>e) ^ Consequently, as 
follows from (jlip and (|2ip . this map satisfies 



Therefore, the map Vs in B{0,3p) fulfills the conditions of Theorem 4.1 in |27]. This theorem 
guarantees that: 

• we may cover the set {x G -6(0, 3p — e/b) \ \vs{x)\ < 1 — (e/ft)^/*^} with a finite collection of disjoint 
balls := {B^j}; 

• the radius of B^, rad(;S'^), which is defined as the sum of the radii of the balls Bj, rad(;B^) : = 
J]^.rad(S|), satisfies rad(i3^) < IQ-^S ■ dist(w,aK); 

• denoting dj = degg^e(f£) if B^ C B{0,3p — e/b) and dj = otherwise; 
we have 



Note that, by the construction of in i?(0,3p) \ B{0,p), if we have degg^e(ue) ^ then i?| C 
B{0,5p/2). Thus dj = degQQi:{ve) for all j. 
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In order to obtain a lower bound for we use the identity 

F,{v,,B{x,,e^l')UB^) = 1/ /|Vz;,|2 + ^(l-K|2)2| (33) 



2 J B{x^,e'^/*)UB^ 
f 



1 

+ 2 



|2 



[{U^-b')\Vv, 

The first integral in ()23p is estimate via (|22p : 

||V..P + ^(1-NY| > vr62^|deg,^^.(.,)|ln^-C 



> TT^^doln- -Co. (24) 



By combining (|20p and Proposition [21 we have for small e 

2 yB(x,,ei/4)ue= I 2e J 

-'-^L. {|V-e|^ + ^(l-klY}-C>(l-6^)|lne|V3_c'; 



(25) 



here we rely on the assumption (jH|) on the behavior of 5{e) as e — )• 0. 

Substituting the bounds (|24p and (|25p in (|23p we obtain a contradiction with (|lip . This com- 
pletes the proof of Proposition \5\ □ 

3.2 Distribution of Energy in _B(0, |) 
Proposition 6. The following estimates hold: 

\l U^,\Vv,\^ = 7^dl\\n5\+0{l), (26) 



and (recall that ^ ~ ^) 

F^(^„i?(0,l)) =7rdo62|lnC|+0(l). (27) 
Proof. We start by proving that 

F5(t)„i?(0,l)) >7rdo6'|lnC| -0(1). (28) 

As before, we use Theorem 4.1 in [2J7j: for < r < tq := 10~^ • dist(ci;, (9i3(0, 1)), there are 
C > and a finite covering by disjoint balls Bf,..., Bf^ (with the sum of radii at most r) of the set 
{x G B{0, 1 - ^b) 1 1 - \ve{x)\ > (e/6)V8} g 



I /^^^{|V^.P + ^(1-|^.P)2} >vrD|lne|-C, 



(29) 



with D = and 

'deggB^iVe) if C B(0, 1 - ^/6) 



otherwise 
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From Proposition |5l for e small, if degQge{ve) 7^ then i?| C -6(0, 1 — tq) C -B(0, 1 — It follows 

that D > do and then (|28p is a direct consequence of (|29p and the bound C/e > 6. 
We next prove that there is C > s.t. 

-/ U^lVvsl"^ >Trdl\ln6\-C. (30) 

2 JB(0,f)\i3(0,l) 



By Proposition |5l {v^l > 1/2 in -6(0, j) \B(0, 1), therefore, := ||^ is well-defined in this domain. 
Observe that 



2 ./b(o,£)\b(o,i) 2 







/ 

/B(0,f)\B(0,l) 


V, „ 



^dgln^. (31) 



We claim that (|30p holds with C = vrdQl lnp| + 1 (for small e). By contradiction, assume (|30p does 
not hold. Then, up to a subsequence, we have 

U U^\Vve\'<^dlln^-l. (32) 

^ JB(0,f )\B(0,1) ^ 



On the other hand, we have 
and therefore 



|V{},|2> / \Vwe\^-l {I - \v,\^)\Vw,\\ (33) 

JB{0,f)\B{0,l) 



Since > ^ in i?(0, |) \ 5(0, 1) we have |Vt&e| < 2|Vi)£|. Therefore, by p2]) . Proposition [5] and 
(iHl) we estimate the last term in 



f (1 - |{),|2)|Vu;,p < Csllneri / |Vi),|2 < Cg^^ ^ 0. (34) 

Js(o,f)\s(o,i) JB(o,f)\i?(o,i) |lne|3 

Combining ([31]), ([331) and (|34|, we find that 

/ |V?)e|^ >vrdgln^-Oe(l). 

JB(0,£)\B(0,1) 



Since \Ue - 1| < C^^ in Br \ B{0, 1) (see Proposition [2|), we obtain a contradiction with (|32p . and 

5 

(|30p follows. Comparing the lower bounds (|28p and (|30p with the upper bound in Proposition H] 
the Proposition [6] follows. □ 

Using exactly the same techniques as in the proof of Proposition [6l one may easily prove the 
following estimate. 

Corollary 1. For any R2 > Ri > 1 



F^{v„B{0,R2)\B{0,R,)) = Oil). 

3.3 Convergence in C°°{K) for a compact K s.t. flU = 
Proposition 7. Let K C M? \Ld be a smooth compact set. Then we have 

Vs is bounded in C^{K) for all k > (35) 

and there is Ck > s.t. 

\ve\>l- Ck^'^ m K. (36) 
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Proof. From Proposition [2] 

E^{U,,K) = \j^ iVf/ep + ^(1 - U'^f = 0{l). (37) 

As in jlH], the following expansion holds 

E^{UeVe,K) = E^{U„K)+F^{v„K) + / {\Ve\^ - l)Ued,Ue. (38) 

JdK 

Using (fTOjl . we have 

/ - l)Ued,lJe = 0.(1). 

JdK 

With (|37p and (j38|) . we conclude that E^iUeVe, K) = 0{1). Since [/. and C/e'Oe satisfy the Ginzburg- 
Landau equation —An = ■^u{l — |np) in K, as well as It/.j < 1 and iJ/^nel < 1. Theorem 1 in |23] 
implies that 

Us and UsVe are bounded in C^'iK) for all k>0. 

It follows that Vi; is bounded in C^{K) for each A; > 0. On the other hand, using the fact that Vs is 
bounded in C^[K) together with the equation of -u., we find that 1 — \ve\'^ ^ CkS,"^ in K. □ 

Corollary 2. For K cR'^\ CO, up to a subsequence, there is some vq G C°°(i^, S^) s.t. — t- vq in 
C^{K). 

We are now in position to bound the potential part of the energy. 
Corollary 3. There exists C > independent of e s.t. 

^/ (l-klY = i/ {l-\ve?f<C. (39) 

Proof. Note that from Propositions |4] and |6l we find that there is C > s.t. 



1 



(l-\vJ^?<C. 



B{0,p/5)\B{0,1) 



Thus it remains to prove the estimate in B{0, 1) for small e. Using (|35p . trQQ(^Q i)V£ is bounded in 
C^{dB{0, 1)) and 1 — \ve\'^ ^ C'C^ on dB{0, 1) (for small e). These properties, allow us to construct 
a smooth extension Vg of trg^(o,i)^£ into 5(0,2) \ 5(0,1), s.t. h = tigB(o^2)^e is S-'^-valued and 
independent of e, 1 - < C^^ ^(g, 2) \ 5(0,1) and 



/ \mf + :^,ii-\vs\r]<co. 

J B{0,2)\B(0,1) I ^? J 



'B(0,2)\B(0,1) 

(For example, this construction is performed by mimicking (|2ip ) 



(40) 



Define We as We = Vs in 5(0, 1) and We = in 5(0, 2) \ 5(0, 1). Clearly, We G 5^(5(0, 2)), 
is bounded in L^(5(0,2)) and, thanks to Proposition |6] and (|40p . 

J / l\Vw,\^ + ^{l-\we\^A <7Tdo\lnC\ + Co. 

^ JB(0,2) I ^4 J 

We may now apply Proposition 0.1 in [14J to in 5(0, 2) to conclude that — 9 / (1 ~ l^eP)^ ^ 
Ci. Therefore the bound (pQl) holds. □ 
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3.4 The bad discs 

Consider a family of discs {B{xi,e^^^))i^j such that (here / depends on e) 

for all i £ I we have Xi G il, 

For ^ G (1/2, 1), let C = C{fi), £o = £o{n) be defined as in the second part of Lemma[T] For e < Eq, 
we say that B(xi,£^^^) is fi-good disc if 

Feive,B{x^,£^/*)nn) < C(^)|lne| 

and B{xi,£^^^) is fi-bad disc if 

Fe{ve,B{xi,£^/^)nn) > C(^)|lne|. (41) 
Let Je = J := {i G / 1 i?(xj, e^/^) is a /u-bad disc}. 

Lemma 2. There is an integer N , which depends only on g and fi, s.t. 

Card J<N. 

Proof. Since each point of fl is covered by at most 16 discs B{xi,£^^^), we have 

Y,Fe{ve,B{x^,e^/^)nn) < 16F,{ve,n). 
iei 

The previous assertion implies that Card J < ^^^^ . □ 
The next result is a straightforward variant of Theorem IV. 1 in 

Lemma 3. Possibly after passing to a subsequence and relabeling I, we may choose J' C J and a 
constant A > 1 (independently of e) s.t. 

J' = {1,...,N'}, N' = Cst, 



>8Xe^/Uori,j eJ',i^^ 



and 

U,6j5(x„e^/^) C Ui^,rB{x^,Xe^/^). 

We will say that, for i G J', B{xi, Xe^^^) are separated fi-bad discs. From now on, we work 

with separated ^-bad discs. Denote Xj = — . By Proposition [5] we know that for small e, we have 


Xi € Bi. Clearly, up to a subsequence, 

there are ai, ...,0^, k distinct points in Bi 
< {Ai, Ak} a partition (in non empty sets) of J' s.t. (42) 
for i ^ J', if i G A^, then xi ^ a^. 

Note that for i G J', we have 

fv?? > 0, for small e, 
2/ G |qi, Qk| <^=^ < (43) 
I there is a /i-bad disc inside B{y,rj). 
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3.5 Convergence in Hl^^i^'^ \ {'^i' Q^k}) 
We have the following theorem. 

Proposition 8. Let ai,...,aK be defined by (j42p . Then we have: 

1. The points ai, belong to uj. 

2. There exists vq £ H^^^{M? \ {ai, ...,Qk},S^) s.t. (possibly after extraction) 

Ve VQ m Hl^^i^ \ {ai, a«}) (44) 



V, ^ m CiO,(M2 \ {a,, ...,aj). (45) 
3. There exists r/o > s.t. for all < r] < ijq and for sufficiently small e we have 

^^^dBr,(ak){Ve/\Ve\) = deggs^^ (^^) (wq) = 1. 

4- K = do. 

Proof Step 1: v^ vq in Hl^iR"^ \ {ai, ...,a^}), vq G /7/„^(IR2 \ {„!,..., a«}, S^) and G w 
Proposition |5] guarantees that ai,...,a«; £uJ. Let 

_ fl0"2 • miufe^fc/ |afc - afc'l if k > 1 
"^""[l ifK = r ^ ' 

Applying Theorem 4.1 in [27] we have for all < i] < r]Q and for small e 

1 f r,_. ,o 52 



I [ \m\' + ^il-\vs\r]>-dolnyC 



(47) 



with C independent of e and ry. Combining (|47p with (|27p and Corollary [T] we obtain that v^ is 
bounded in H^[K); here IT C \ {ai, is an arbitrary compact set. Therefore, there exists 
Vq G -^[^^(M^ \ {ai, ...jOk}) s.t. we have Ve vq in Hl^^{E? \ {ai, ...,0^}) (possibly passing to a 
subsequence). Since ||1 — |we|||L2(-^) — t- for all compact sets K CM.'^ \ {cti, ...,0^}) we find that vq 
is S"^- valued. 

Following the proof of Step 7 in Theorem C in |18j, we can prove that ai,...,an ^ duj, thus 
«!, G UJ, and the first assertion follows. 

Step 2: Proof of 2. 

Adapting the techniques of [lOj (Theorem 2, Step 1), we establish (|44p and (|45p in a ball 
B = B{y, Ro) s.t. BcR^\{ai, a J. 

Let y G ]R2 ^^^^ let R' > R > be s.t. B{y,R') C \ {a^,. Since F^{vs, B{y, R')) 
is bounded independently on e, there is Rq G {R,R') (independent of e) s.t., passing to a further 
subsequence if necessary we have 

/ j \drVe\'^ + 72 (1 - 1*^1^)^ \ ^itl^ independent of e. (48) 

JaB{y,RQ) I C J 

Indeed, for r G {R,R') denote 

Is{r)= I \\yie? + ^,{l-\Ve\'f 



dB(y,r) 
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Using the Fubini theorem and the Fatou Lemma we have 

f-R' f-B' 



0< / hminf /e(r) dr < Uminf / Is(r)dr<C'. 
Jr ^ ^ Jr 



C 

Consequently, hminf£/£(r) < oo for ahnost all r G {R,R'), so that (|48p holds with C 



R' -R 

Let = tigsVe- Since \vs\ > 1/2 in B = B{y,RQ), we have degQQ{g^) = 0. The bound (|48p 
implies that, up to choose a subsequence, Qir is weakly convergent in H^{dB). Consequently there 
is he H^{dB,S^), h = e"^, ip G H^{dB,M.) s.t. 

Qe ^ h uniformly on dB^ (49) 

Qe^hm H^''^{dB). (50) 

Let r^e : S — 7- M"*" be the minimizer of / <^ |V?7|^ + -^{l - r]f \ in H} |(S,M). Then ry^ satisfies 

+ = 1 in 5 
\% = \9e\ on dB 

It follows from |1U) that 

/^{|V7?e|' + ^(l-%)'}<Ce. (51) 

Using (|49p . there is (^g G H^{dB s.t. = l^ele**^^ and (/?£—)• 93 uniformly on (95. Following 
|10) . denote by G -?/^^(i3,M) the unique solution of — div(a^V'0e) = 0. (Here a = 6 in a; and 
a = 1 in \ lj.) From (|50p . — )• ^/^ in H^{B) where tp G H^{B,M.) is the unique solution of 
-div(a^VV') = 0. Since rjee''^^ G Hg^{B), we have 

F^{v„B)<F^{7j,e'^^,B) < \ j C/^IV^.^ + Ce ^ J / a'\VM\ (52) 

On the other hand, since vq in H^{B), we have tig = e*"^ with (/> G H^{B,M.) and 

liminfF5(i)e,5) > liminf J / ^/^iVw^p > ^ / a^lVvoP = J / a2|V(/>p. (53) 
= ^ ^ Jb ^ Jb ^ Jb 

(The last inequality follows from [/^ — )• a in L^, \Ue\ < 1 and Vs vo in -fT .) 

By combining (|52p . (|53p and the fact that ip minimizes / a^|V • |^ in H^{B,M.), we find that 

(|44p holds. Furthermore, the map ^ in (|52p is the same as (j) in (|53p . 
Note that since ^ 

^ ^ "''^ o,(l)<F^(t)„i?), 







IB 


l^el 



2 

by comparing (|52p with (|53p . we also have 



V\Vef + hl-\Ve\'f^^. (54) 
K ? 



In order to prove (|45p . it suffices to establish the convergence 

^ in L°°(S) with cpe G i/^^_(5,M) and = \ve\e'^% (55) 
and to use the fact that — )• 1 uniformly. 
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Proof of (|55p . If du} n S = 0, then the argument is the same as in [lOj. Assume next that 
3a; n i? / 0, and let tjj £ H'^^'^{B,R) be the harmonic extension of ip. Since C — V' S Hq{B,R) 
satisfies — div(a^VC) = div(a^V'0), Theorem 1 in [22] imphes that (f) G W^''p{B,M.) for some p > 2. 

We next prove that, for some q > 2 and B = B{y',R) s.t. B{y',2R) C B, we have — 
^llvFi''J(i?) ~^ (Oiice proved, this assertion will imply, via Sobolev embedding that (|55p holds.) 

Note that (up to a subsequence) i;^^ — )• in L^(i?,]R). Thus we have 



-div 



div 



(C/||i),|2-a2)V0 



in 5 



V2|~, 1 2 

V 

From Theorem 2 in [22], there is 2 < g < p and C > s.t. 

l|V(c/>. - C/.)||^,(5) < C (i?-2+2/5||^^ _ ^11^^^^^ + ||(f>2|^^|2 _ a2)V</.||i,(5) 

Consequently, \\4>e - 4>\\wi,g(^B) 0- 

Step 3: We prove the third assertion 

Let rjQ > r] > 0, with r/o defined by (|46p . Denote dk = '^^SdB{af,,r)i'^o)- These integers do not 
depend on r G {ri,r]o). Moreover, we have Ylk^k ~ '^0- For r G {r],rjQ), we obtain that 



0. 



27r|4| < 



dB{ak,r) 



1/2 



|(?rfo| 



and therefore 



iVvor > vr^ln^. 

BK-,77o)\BK,r?) ^ 



Consequently, we have 



lim inf 



VJkB{ak,m)\B{ak,ri) 



IV7- |2 ^ ^ 
I £1 - 2 



> vr In 



UfcB{Qfc,»y())\B(afc,7j) 



|Vt;or 



(56) 



By combining (|47p and (|56p . we obtain the existence of C independent of e and r/ s.t. 



Ufc-B(afc,?;o) 



> vrVdlln-^ + ^dolnT-C 



Trdo 111 + ■^(^^ ^fc ~ '^o) In C. 



Therefore, dk must be either or 1. Otherwise, (|27p cannot hold for small r/. Applying the strong 
convergence result from Step 2 with K = B{ak,rf) \ B{ak, 5), we have that for small e. 



dk = degg5(Q.^ 



We next prove that = 1 for each k. By contradiction, assume that there is ko s.t. dk^ = 0. We 
may assume that ko = 1. From (|43p . there is a (separated) |U-bad disc B(xo, Xe^^^ /6) in B{ai,r]o). 
Thus by (|4ip . we have 

F^ive,B{xo,Xe'/y5)) > Cif,)\lne\. 
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On the other hand, since in \vs\ > 1/2 in B{ak,riQ) \ i?(afc, ?7o/2), applying Theorem 4.1 in |27) in 
B{ak, t]q), k G {2, «;}, with r = 10~^ • rjo we find that 

F^{ve,B{ak,m)) > &Vegai?K,,,o)(^e)l|ln^| - C, k = 2,...,K. 

Since X]fc=2 '^^S9S(«fc,»7o)(^e) ~ ^O' above estimates yield 

F^{v„B{0, ^)) > b^do\ In^l + C7(//)| lne| - C 

which is in contradiction with ([H]) and (|lip . Thus dfc = 1 for k £ {1,...,k;} and consequently, 
K = do. □ 

We are now in position to estimate the rate of uniform convergence of l-u^l in a compact set 
K cR^\{ai,...,ad,}. 

Corollary 4. There is C > s.t. for r]Q > ij > and small e we have 

\ve\ > l-C|lne|"i/3 in B(0,^)\B(ai,7j). 





Proof. Due to (|36p . it is sufficient to establish this result in -6(0, 1) \ B{ai,r]). Combining Corollary 
[3] with (|44p . we obtain that 



F^{v„B{0,2)\B{a,,ri/2)) <C{ri). 
Thus for all x G B{0, p) s.t. B(x, e'^/'^/6) C B{0, 2) \ B{ai,r]/2), for small e we have 



F,{v„B{x,e^/^)) < Fi:{v„B{0,2)\B{a„v/2)) < | In 



1/3 



From Lemma [T] (first assertion), we obtain the existence of C > (independent of e and rj) s.t. 
|f£(x)| = |v£(x)| > 1 — C| lne|~-^/^. Finally, since for all x G -6(0, l)\B{ai,ri) we have B{x, C 
B{0, 2) \ B{ai,r]/2), Corollary 1 follows. □ 

3.6 Information about the limit vq 

Following |11) (Appendix IV, page 152) we have 
Proposition 9. For all 1 < p < 2 and for any compact K C M^, is hounded in W'^'''^{K). 

Let 9i be the main argument of |^~°'| and set = + ... + O^q. Note that V0 is smooth 
away from {ai,...,aK} and njj|5^ = e*^- Let g := trgBiVo and cpo G C'^{dBi,M.) be s.t. g = 
Hi ifEf^e^'^o = e*(^+'^o) (see |l2] for the existence of tpo). 

Proposition 10. The limit vq satisfies — div (a^ f o ^ Vf^o) = m T>'{M?). Moreover we may write 
Vq = Here 93^ is the solution of 

{ -div [a2V(0 + if,)] =0 inBr ^^^^ 
= V^o on dBi 

Proof. Let (/> G P(M^), and set K = supp((/)). By Proposition El we have {7| x ^ VtiQ 



in LP{K) for p < 2. Multiplying the equation —div 
parts, we obtain 



by (/) and integrating by 







/ —div 


Ve X VVe 


Ik 





I Ve X VVe ■ V(t) 

Jk 



a Vq X Vvq • V0 = / —div (a vq x V^o) cf). 
K Jk 
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Consequently — div [a? vq x Vfo) = in P'(M^). 

In order to prove that -div [a?V{d + tp^,)] = in \ {ai, a^J, we follow [H], Step 12 of 
Theorem C. 

Next, we prove that ip^, is harmonic in a neighborhood of a^. Fix A > and xq ^ oo s.t. 
B{xq,2\) d 00 \ {ai, Ofig}. As we estabhshed in Proposition [8l Step 2, F^{ve,B{xQ,2\)) is 
uniformly bounded in e. Proceeding as in Step 2, we conclude that exists Aq G (A, 2A) s.t., after 
passing to a further subsequence, we have 

/ ||V^.P + 4(1-N')'|<^ (58) 

JdB{xoM) I J 
with C and Aq independent of e. Now, if minimizes 

Edu) = U ||Vn|^ + :^(a^-NY 
2 yB(o,f ) I 2^ 

subject to -6(2;) = fs{Sx) on dB(0, 0, then minimizes E^{u, B{xq, Xq)) with respect to its own 
boundary conditions. In other words, := ^ minimizes the classical energy 

II ||V^e|^ + ^(l-|-e|Y 

2Jb(xo,Ao) I 2^ 

among w G H^{B{xo, Xq)) such that w = /i^ := ^ on dB{xo, Xq). It follows from (|58p and 
Proposition [2] that /i^ also satisfies 

/ {\drK\' + ^{l-\h,\'A<C + l. (59) 

JdB{xoM) I 24 J 

Note that by Proposition [2] we have 

\\we\\Loo(B{x,M)) < l + ce-'^o^. (60) 

Using (|60p and the uniform bound from Corollary IH we may repeat the arguments of Theorem 2 
in |10) and conclude that, up to a subsequence, there exists an S^-valued map wq s.t. for every 
compact K d {ai, ...a^^p}) we have 

We wo in C^iK) (61) 

and 

^ |V^i;o|2 in C-(K). (62) 

„. . fminlafc — a,| d[st(ak,^io)^ , , , t / o n t 

l^ix now r < mm < > and denote Ur := |x G uj,dist[x,ouj) > r\. It 



follows from (jbip that Ws — ?• qo '■= trg^^wg in C°°{duJr)- In view of Proposition O we have wq G 
VK^'^(wj.), p <2. By Remark I.l in |11] . this implies that 

Wq = w exp Cfc In I x — a/j I + ix^ . 

Here: 

• z/; is the canonical harmonic map (see |11] . Sec. 1.3.) having singularities {a^, A; = 1, ...,do} and 
equal to qq on So;,.; 
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• the Cfc's are real coefficients; 



• X is the solution of 

Ax = in ojr 

x{x) + ^k'^k^'n\x - ak\ = on du>r 

Repeating the argument of |11) . Theorem VII. 1, Step 2 (the key ingredients of this proof are (I6ip . 
(|62p and Corollary |3|) , we find that = 0, A; = 1, do, and, consequently, wq = w in cOr- Finally, 
by [llj . Corollary 1.2. , we know that the canonical harmonic map w is of the form w = e*^^"^'^*-* with 
(p^ harmonic in ujr- 

□ 

3.7 Uniqueness of zeros 

Proposition 11. Fore sufficiently small, the minimizer has exactly do zeros. 

Proof. It suffices to prove that for small e there is a unique zero of in B{ak, r), k = 1, do, with 
r defined in the proof of Proposition [TOl 

Since Ws = from Proposition |2] and Proposition [10] we see that wq = vq = e^^^'''^^''^ in 

B{ak,r), where 9k is the phase of ||^~"^ and = (/j^, + -0^ is harmonic in B{ak,r). Using (|6ip 
and (|62p and arguing as in the alternative proof of Theorem VII. 4 in | TT| (page 74) we obtain that 
VHkiak) = 0. 

Finally, we are now in position to obtain, as in Theorem IX. 1 |llj (using the main result of IlTJ), 
that there is a unique zero of Ws (and, therefore, of Vs) in B{ak,r). 

□ 

3.8 Summary 

We have thus proved 

Theorem 3. Let Sn i and Ve„ be a minimizer of pzp in (jlSp for e = e„. Then there exist do 
distinct points ai, ...,0(1^ G and a function vo ^ Hl^^{M.'^\{ai, ...,aci^},S^)r\W^^^{M.'^,S^) (p < 2) 
s.t., up to a subsequence 

1. Vsn^vo in H^^^iR'^ \{ai, ...,ado}) and Ci° \ {ai, arfj), 

2. Ve„ - vo m W^^^{R^) (p < 2), 

3. forK <^R^\{ai,...,ad,}, jw^J > 1-| lne„|-i/3 m K and [ \V\vsJ\^ + ^{l-\ve,ff^O, 

Jk ? 

4. for K (s R2 \ ve„ vo m C°^{K) and 1 - \vej < Cr^'^, 

5. Ve„ has exactly do zeros x", ■■■,x2^ and — )• at, 

6. Vo satisfies — div {a? vo x Vvo) = in D'(M^). 
Let us summarize the proof of Theorem |3J 

• Statement 1. is established in Proposition |8l 

• Statement 2. follows from Propositions |9] and [101 

• Statement 3. is a consequence of Corollary [4] and (|54p , 

• Statement 4- is Corollary [21 
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• Statement 5. is proved in Proposition 1111 

• Statement 6. is established in Proposition [101 
The proof of Theorem [3] is complete. 

4 Renormalized energy for the model problem 

In this section, we establish the expansion for Fi;{v£, B{0, p)) = F^{v£, B(0, |)); specifically, we 
derive the expression for 

lim^F^{ve,B{0,^)) -Trdl\ln6\ -7rdo6^|lnC|} • (63) 

(We are going to prove that this limit exists). 

In order to find an expression for (|63p . our strategy is the following: 

• (Section 14. ip We first study the minimization of the Dirichlet functional among S^-valued maps in 
annulars B{0, p/5)\B{0, 1) with the Dirichlet boundary conditions: f^{S-) on dB{0, p/S) and on 
dB(0, 1). Here = where is given in the model problem and g^,g^ £ C°°{dBi,§^) are s.t. 

g^ — ;> (7° in C^. We get that the Dirichlet energy has the form Trdg ln{p/6)+WQ{fo)+Wi{g'^)+os{l). 

• (Sections 14. 2| |4?3] and l44l) In B{0,1), we study the weighted Ginzburg-Landau functional with 
the Dirichlet boundary condition g^ on dB{0,l). Making use of the previous bullet point, one 
may obtain the matching upper and lower bounds and use them to derive the third term of 
renormalized energy, which depends on the limiting locations of the zeros /3 = (/3i, ...,/3rfg) E uj'^° 
and on g^. We establish that 

inf F^{; B{0, 1)) = ndob^ In ^ + dob^j + miP, 5°) + o,(l). 

a" 

• (Section 14. 5p Finally, we make a fundamental observation: the limiting function go = limtrg^jUg 
and the points a obtained form Theorem |3] form a minimal configuration for Wi{g) + W2{(3,g). 
Thus, introducing 

W{p)= inf {WiCg) + W2{P,~9)} 

with deggg^ {g)=do 

we conclude that a minimizes W. 
In this section we prove the following theorem. 
Theorem 4. The following energy expansion holds when e — t- 

F^{ve,B{0, ^)) = TTdob^ In ^ + TTdl In ^ + VFo(/o) + W{a) + d^b^^ + o,(l). (64) 

Here the points a = (ai, ...,ado) are obtained from Theorem\^ 'j > is an absolute constant and 
WQ{fo),W{a) are renormalized energies: 

• Wq is independent of the points ai, and given by < I7^) . 

• W is given by (|9Up . it is independent of fQ and the limiting points (ai,..,ad(,) minimize W. 

Remark 2. The renormalized energy in the expansion (|64p decouples into the part that depends 
only on the external boundary conditions Wo(/o) and the part that depends only on the location of 
the vortices W{cx). Since a minimizes W, the external boundary data has no effect on the location 
of vortices inside the inclusion. This is a drastic difference with the results of [11] and [18], where 
the Dirichlet boundary data on the external boundary influences the location of the vortices. 
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4.1 Minimization among S^- valued maps away from the inclusion 

Denote Bp := B{0,p). Let 

. 1 n 1 if^ fmC\dBp) 

{f%^s<iCC°^idBp,S'), f eC^{dBp,S')hes.t. \ S^o 



and 



.deggBi(5 ) = do 



For 5 £ (0, 1), we denote As = Bpj^ \ Bi and 



W5 = {u(^ H\As,S^) I troB^^X-) = f^i^-) and troB.u = /}, 

Ys = {uG H\As,S^)\tT9B^^X-) = fi^-) and trQB.u = 5°}. 
Consider the following minimization problems: 



Is{f',g') = Is= inf I [ \V 

Js{f,g') = J5= inf I [ \V 
u£Ys 2 J A, 



Proposition 12. For small e, Is is close to Js, namely 

Is = Js + 05(1). 



(Ps) 

m 

(65) 



Proof. In this subsection 6 stands for the main argument of z ie. ^ = ^ . For 5 > 0, let 
^S G C°°(95i,M) be s.t. / = e'^'^^e+M and Q G C°°{dBp,^) be s.t. /-^ = e<'^<i^+^s) ^ We may 
assume that </)5 — )• 0o in C^{dBi) and (^5 — )• ^0 in C^{dBp). Note that 

« G ^ -u = e<'P+'^^^^) with 99 G (66) 

Here ws := W ^ H^{As,M) \ iTQBpf{-) = C<5(5-) and iiaB^V = 

s 

Since = in ^5 and d^O = on dAs, for n G Ws we have 

/ |Vn|2= /" \v{ip + doe)\^ = 4 [ \ve\^+ [ \Vip\^ 

J As J As J As J As 

Consequently, the problem (Ps) has a unique solution us = e'^i'^od+'Ps) ^ with (/p^ being the unique 
solution of 

— Aips = in ^5 
< ¥'5(-) = Cs (5-) on 5B£ . 

o 

^ = <?!)5 on (9^1 

With the same argument, the problem (Qs) admits a unique solution vs = e^i'^od+i^s) with tps being 
the unique solution of 

'-A^S = in As 
M-) = Co{S-) ondBE. 

o 

,"05 = 00 ondBi 
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Denote rjs = tps — ips- Then ijs is the unique solution of 



Ar]s = in ^5 

% = 6 - Co on OBe . 

o 

= fps-fpo on dBi 
(Here C{x) := C{6x)). 

One may prove that ||^5||l2(^^) is bounded and more precisely we have the following result. 
Proposition 13. 



As 



|</'o|^i/2(§i) + ICo|^i/2(9B^), as 



5^0. 



Proof. Let (a„)„gz, (&n)nGZ C C be s.t. 



m9 



We have 



(e^^) = and (oipe'') = ^ 6ne 



From [8] (Appendix D.), denoting -R((5) = |, we have 



1 
2^ 



n 



lni?(5) ■ i?(<5)2M _ 1 



|a„P + |6nP)(i?(<5)2H + i) 

-2(o;76„ + aX)i?(5)l"l 



I (/'O I ^1/2(9^1) + ICoIhi/2(OBp) + 05(1). 
Consequently, as 5 — t- 0, we obtain (|67p . 

Following the same lines as Proposition [13] we obtain 

II V995||j;^2(^^) < C with C independent of 5, 

and 



(|a„P + |6„|2)-(a^6„ + a„5„)/?(5)H 



l|V%||i2(^^) ^0. 



It follows from 1^ and (|69| that 
72 



^ J A., ^ J Ax J A, ^ J As 



I As 

Js + osil). 



(67) 



□ 

(68) 
(69) 



(70) 
□ 
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From (1701) and (1671) . we deduce that 



-'(5 — -"5 -I- '-'SK^J — "t"'^ ^ 

with 



/5 = J5 + 0^(1) = ndl In I + H^o(/°) + + (71) 



Woih) = KolW^dB,) and = ml./2^sB,y (72) 

One of the main ingredients in the study of the renormahzed energy is that the Dirichlet condition 
/min(a:;) = TOjfp^, 70 e minimizes Wq. More precisely, for all /o G C^(5Si,S^) s.t. degg5^(/o) = 
do, we have 

VFo(/mm) = < Woifo). (73) 

4.2 Energy estimates for S^-valued maps around the inclusion 

Let £ C°°{dBi,§^) be s.t. degg^^(g'^) = do > 0, /3i, ■.■,f3(ig are do distinct points of to, 

rjo := 7min|dist(/3i,9u;),min|/3j -/3,| 

For r G (0, ?7o)i we define 



a, :=5i\UfcS(/3fc,r), 
:= {li G i/^(Or,§^) I trasi^ = 5° and dega5(^^ ,,)(n) = 1} 

and 

X - f3i 



Tr '■= G H^{i}r, S^) I trgg^ti = and there are ■ji G §^ s.t. tr^^^^, ,,)n(x) = 7i- 



|x - /3i 

Consider two minimization problems 



K{r,g\l3) = K(r) = inf \ f a^\Vu\^ (Rr) 



and 



Lir,g\P) = L{r)= ini \ [ a'\Vu\\ (3 = {Pi, Pd,] ■ 



We denote 9 = 9i + ... + 9^^ where 9i is the main argument of -j — , i.e., 



\x Pi\ \x Pi 

Let ■00 be the unique (up to an additive constant in 27rZ) solution of 



-div [a2(VV'o + V6I)] = in Bi 
g»(e+^o) = ^0 Qj-^ 



Lemma 4. ([18], Appendix A.) 



with 



(Sr) 



(74) 



K{r) = ^j a^\V9 + Vi}o\'^ + 0{r\ In r\), 
Hr) = ^j a2|V^ + VV'o|^ + 0(r|lnr|), 

If a^\V9 + V'iPo\'^ = 7Tdob'^\lnr\+W2i(3,g^) + 0{r'^). (75) 

In (|75p . iy2(/3;5'^)) whose explicit expression is given in |18] . formula (106), depends only on /3 
and g^. 



22 



4.3 Upper bound for the energy 

Lemma 5. Fix p > and let G C°^{dBp), fo G C~((9Sp,§i) be s.t. fo in C\dBp). Let 

p = (/3i,...,/?do) e w'^" be s.t. Pi 7^ Pj for i ^ j. Then, for each G C~(5Bi,§^), the following 
upper bound holds: 

inf F,<ndobHn^ + 7rdlln^ + Wo{fo) + Wi{g'') + W2{P,g'') + dob^-/ + o,{l). (76) 
Here Wq, W\ are defined by ([72|) and W2 by (j75]) . 

Proof. We construct a test function G Hj {Bpj^, C) which gives (|76p . Fix < r < ryo- Let 



Us be the minimizer of [Ps) with (7*^ = and f^ = rj 

III 



and 



be the minimizer of (1^ 



Note that /"^ — )• /o = hm^ in C^{dBp). For each i = l,...,do let be the global minimizer of 
the classic Ginzburg-Landau energy in B(f3i,r) with the parameter and the boundary condition 
u^'^{x) = h^{x) := on dB{f3i,r), 7^ G is defined through Ur- Denote 



I(e/6,r) 



inf — 



i/ { 



vnfn^ + ^a-K'T) 



Vn|^ + ^(1-HY 



(77) 



Lemma IX. 1 in [llj implies that 



6r 

/(e/6,r) = ^ln- + 7 + 05(l). 



(78) 



We next extend the tij's to -Bp/5. For this purpose, we consider G C°°(M, [0, 1]) s.t. ^ = in 
and = 1 in [1, 00) and set 



X.(se*') = C(s-7 + l) \fe\{pn-^ 



+ 1. 



In view of (|14p . we have \\xe — M\L'^(Bp/g) — C'e. Consider the following test function 

'XeUs mBp/s\Bi 



u = < 



Ur 



in Bi\UB{(3i,r) 



(79) 



^uj' in B{/3i,r) 



Clearly, 



inf F, < F^iwe) < ndl In ^ + ^^1(5") + W2iP, 5") + W^o(/o) + ^rdob'' In ^ + dob''j + o,(l) + /i(r) 
with h{r) = 0r{l). Thus, letting r — )• as e — )• we obtain the desired upper bound. □ 
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4.4 Lower bound 

We prove that the upper bound (|76p is sharp by constructing the matching lower bound. 

Lemma 6. Let e„ \, 0, v^^ be a minimizer of (|17p in (|18p for e = e„ and a = (ai, ...,adQ) G 

be given by Theorem\^ Denote g^ := hmtrg^^^e G C°°((?i?i, S^). Then the following lower bound 

holds: 

F{v„Be) > vrdo6' In ^ + Trd^ In ^ + Woifo) + Wi{go) + W2{a, 9o) + do^S + Oe{l). (80) 



Proof. As in the proof of Lemma |5l we spht Bp into three parts: Bp \ Bi, Bi\ UB{ai,r) and 

5 6 

[jB{ai,r) with small < r < r/o- 

In i?£ \ U-B(aj,r) one may write = \ve\we. Using Corollary |4] and (|44p we have 

6 



'Bp\-Bi 

> \ l _U^eWWe? + Oe{l). (81) 

1 



We take 17'' = - — and f^ = -j — Note that with this choice of f^,g^ one may apply the 
results of Sections 14.11 and 14.21 From (|8ip we obtain the lower bound in Bp \ Bi : 

S 

F^ive,BE\B;)>J5 + Oe{l) (82) 
with Js the energy associate to the minimization problem i\Qsl (see page I20p . 



Let vq be defined by (|44p . Since we have — >• vq in H^{Bi \ UB{ai,r)) and t/g — )• a in 
L?'{Bi \ UB{ai,r)), from Proposition 1101 and Lemma |5]we obtain 



F^{ve,Bi\UB{a,,r)) >\( a^\\Ivo\'' + o,{l) 

^ JBi\UB{a,,r) 

> If a^\Ve + V^Po\^ + Oe{l) 

^ JBAuBiaur) 



Bi\UB(ai,r) 

= K{r) + 0{r\lnr\) + Oe{l), (83) 

where K{r) is defined by <\Rr\l (see page 122 p . 

In order to complete the proof of the lemma, we need to obtain a sharp lower bound in each 
ball B{ai,r). Actually we will prove that 

Fi:{v„ B{ai,r)) > b^Ii^b, r) + o,(l) + 0^(1), (84) 

with I{S,/b,r) being defined in (|77p . The estimate (|84p is equivalent to 

F^{v„B{ai,r)) > b^I((/b,r + r^) + Or{l) + o,{l). (85) 

Indeed by (|78|) we have /(^, r + r^) - /(^, r) = 0^.(1). 

We now make use of the construction by Lefter and Radulescu in [20] and |19) . From Proposition 
[TOl we know that vq = e*^^'"'"'^*"'"'^') with (^^, ipi harmonic, and therefore smooth in B{ai,r]) {rj > r 
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small). Set fij = ■ Without loss of generality, we can assume that a, = and cTi(O) 

Consequently, |crj(x)| < C\x\ with C independent of r] and < r]. Let 



= X^e<^^+''l) where 
From Proposition |8] and (|54p , we obtain that 

al (Tj in i?^(i?r+r2 \Br), 



Xe 1 in H\B,.^,.2 \ Br) and / _ (1 - ^ 0. 



(86) 
(87) 



Let 



Veise 



1- A. 



(s - r) + A, 



if s G [0, r) 

. g _|_ _|_ \ ^ 

exp <; ^ ( + (T* 5 ]) if s G [r, r + r^] 



Clearly, fi^ G (-B,.+,.2). Consequently, 

62/(^/6, r + r^)< F^iv^, Br) + F^{/3„ \ B~r) + o,(l). 

From (|87p . we easily obtain that 



I _||V|/3,||2 + 1(1 -1/3,1)4=0, 



(!)• 



It remains to estimate 



From (l86 



/ -^{ 



'B^+^2\Br 



-s + r +r 



-s + r + r 



'B,. + ^2\Br 



-s + r + r 



+ Os{l). 



Since |(Tj(se*^)| < Cs, \ds(Ti\ < C and |50-CJi| < Cs we have 



-s + r + r 



cr,- 



-s + r + r (Ti 



< C7 [(1 + r-2) + =0(r 



-2\ 



1 + -2 



Since \Br^r'^ \ = 0{r^) we find that 



Vl9,+ a] 



—s + r + r 



0{r) 



It follows that F^(P£, Brj^r'^ \ Br) = 0{r) + Oe(l). Consequently, (|85p holds and thus we obtain 
Combining (|82|), (|83|) and ([84)), together with dm and we obtain 

F5({}„i?£) > Is + K{r) + b^I{ab,r) + o,{l) + Or{l) 

= TTdl In ^ + 7rdo6' In ^ + W^o(/o) + W{a, go) + 
t. 



+ dob'^-f + Oe{l)+Oril). 



The conclusion of the Lemma follows by letting r — t- as e — )• 0. 



□ 
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4.5 The function qq and the points {ai, ado} minimize the renormahzed energy 

In the previous section, we obtained an expansion for the energy Ft{v^, Bp ) of the model problem. 
To summarize, using (|76p . (|80p and Theorem [3] we get that there are qq = limtrgsj^e and a = 
(ai,...,arfo) e w'^o s 

F^a, Be) = TTdob^ In ^ + vrdg In ^ + W{a, go) + T^o(/o) + ^o^'t + (89) 
s " 

with 

Ty(Q!,5ro) = ^^1(50) + W'2(Q:,fi'o)• 
The goal of this section is to underline an important property of the points a, namely, that they 
minimize the quantity inf^ogj^oo^-g^^ gi) W{-,g^). 
We have the following 

Proposition 14. Let = P^o) £ w'*" be a do -tuple of distinct points and let g'^ G C°°{dBi,S^) 

be s.t. deggQ_^{g^) = do. Then 

W{a,go)<W{f3,g''). 

Proof. Let {P,g^) be as in Proposition [14] Using the test function given by (|79p . we obtain that 
for all e > and r > (small) there is We G {Bp,C) s.t. 

F^iwe) = TTdob^ In ^ + TTdl In ^ + W{P, 9°) + Wo{fo) + rfo^S + h] + h^, 

here h] = 0^(1) and hj = 0{r). 

On the other hand, taking into account the minimality of and (|89p we have 



Wi(3,g'^)>Wia,go) + o,{l) + hl 

The previous estimate implies (letting e — t- and r — )• 0) that W{P, g^) > M^(a, 50) which completes 
the proof. □ 

Thus, for /3 = (/3i, ...,/3rfJ G a;'^^ we define 

W{f3)= inf Ty(/3,g)= inf ^^1(5) + W^2(/3, 5) (90) 

3Si (9)='^o degga^ {5)=rfo 



with M^i and W2 given by (|72p and (|75p respectively. It follows that for a given by Theorem [3] and 
go = troBiVo: 

Wia) = W{cx,go) < WiP) for all (3 = {f3i, ...,pdo) G cu^"- 

5 Proofs of Theorems [1] and [2] 

In this section Ve is a minimizer of in Hg{Q,C). We split the proofs of Theorem [1] and |2] in 
three steps: 

• (Section lS.ip Using estimates on , we first localize the vorticity to the neighborhoods of selected 
inclusions. Then we find two separate energy expansions in two sub-domains of Q: away from 
the selected inclusions and around them. 

• (Section 15. 2p We study the asymptotic behavior of V£. We prove that, for small e, has exactly 
d zeros of degree 1. 

• (Section 15. 3p We give an expansion of F^{ve) up to 0^(1) terms and relate the choice of the 
inclusions with vortices to the renormalized energy of Bethuel, Brezis and Helein. 
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5.1 Locating bad inclusions 

The following result gives a uniform bound on the modulus of minimizers away from the inclu- 
sions. 

Lemma 7. There exists C > s.t. for small e we have 

1. \ve\ >l-C\ lne|-i/3 inn\ ufi^B{ai,6), 

2. there are at most d points Ui-^, ...,ai^, (1 < d! = d'^ < d) s.t. {jfel < 1 — C|lne|^^/'^} C 

Proof. Using Lemma [1] with x = Uii^l"^^^) we obtain that there exist C, Ci > s.t. for e > 
small, 

a Fe{ve,B{x,e^/*)) < \ lne\^ - Ci then \ve\ > I - Cx mB{x,e^/^). 



We prove 1. by contradiction. Assume that, up to a subsequence, there is G \ ^f£iB{ai, 5), 
s.t. |t>e(xe)| < 1 — C| lne|^^/^ with C given by Lemma[T] From Lemma [T]and Proposition [2] 



I [ ||V.,p + -L(l-|,,|Y|>|lne|V3_o(l). 



(91) 



Fix a bounded, simply connected domain Q' such that C 0,', and extend by a fixed smooth 
S^-valued map v in $7' \ Q, s.t. v = g on d^l. 

In view of (1111) for Case I or (1121) for Case II, there exists C > s.t. for small £ 



^/jV„.|^ + ^,I-M^)'<<?|.n.|. 



Therefore, the map Vg in Q' satisfies the condition of Theorem 4.1 |27) . This theorem guarantees 
that 

• there exists = {Bj}, a finite disjoint covering of the set 

{xen'l d[st{x,dQ') > e/b and \ve{x)\ < 1 - {e/bf/^}, 

• such that rad(^^) := rad(5|) < 10^^ . dist(w, dB{0, 1)) • 5, 

• and, denoting dj = \degQQ.{v£)\ if i?| C {dist(x, > e/b} and dj = otherwise, we have 

= Tr^d^llnJI-C, (92) 

j 

with C independent of e. 

Note that since |ue| = 1 in fi' \ 0, if dj ^ then Bj C {dist(x, > e/b}. Consequently, we have 
dj = \degQse{ve)\. 

Assertion 1. follows as in the proof of Proposition [S] (use (|9ip . (|92p instead of (|2Up and (|22p ). 
The proof of Assertion 2. of Lemma [7] goes along the same lines. □ 

We next obtain the following lower bounds for the energy. 
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Lemma 8. For k G {1, ...,d'}, we denote dk = d% = deg^^j-^, ^^-((fg). There exist C,r]Q > s.t. for 
small e and p £ [26, r/o] we have 



2 

and 



If \Vv,f>7TY,dl\lnp\-C (93) 

^ Jn\U? ,B(a,, .o) , -, 



Fe{vs,B{ai„26))>7r\dk\b^\ln^\-C. (94) 

Proof. Let 7?o = 10^^ miiij {dist(aj, 9^2), miiij^j |aj — aj|} and < p < tjq. 

We prove (|93p . By Lemma O > 1/2 in Q\U'I.^^B{ai^,, p), therefore, Ws = is well-defined 

in this domain. From direct computations in 5(ajj,,r/o) \ B[ai^,p) we have 



\l |V«;,p>vrj;d,nn^. (95) 

to a subsequence we have: 



We claim that the bound (|93p holds with C = | lnr/o| + 1- Argue by contradiction: assume that up 



^ ^ 'V^;,|2<^^dfln^-L (96) 



d 

On the other hand, we have 

\VVe\^ = \Ve\'^\VWe\'^ + \^\Ve\\'^ 

and therefore 



f l^Vel^yf IVWel"^ -{l-\Ve\'^)\VWs\'^. (97) 



Using the fact that > ^ in $7 \ uf^^B (oi^, , p) we see that iVf^gl < 2|Vfe|. Therefore, by ([96 
([H]) and Lemma [7] we estimate the last term in ([9 



[ {l-\ve\^)\Vwe\^ <C\lne\-'^ [ \Vve\^ <C^^^ ^0. (98) 

Jn\u'l:^^B{ai^,p) Jn |lne|3 

By combining (|95p . (|97p and (|98p . we see that (|96p cannot hold for small e; this implies (|93p . 
We now prove (|94p . Performing the rescaling x = ^''^ , we obtain 

F,(^;,i?(a,,,25)) = F^{v,B{0,2)) = \ f {u^AVv\^ + ^U^{1 - \v?f] dx, 

where, as in the model problem we set v{x) = v{5x) and ^ = ^. 



By Theorem 4.1 |27) . for r = 10~^ there are C > and a finite covering by disjoint balls 
Bi, Bn (with the sum of radii at most r) of {x £ B{0, 2 - i/h) \ 1 - \ve{x)\ > (^/6)^/^} and 



2 

Dk = Yjj and 



(99) 



degas^. (t)^) if dist(B,-, 95(0, 2)) > e/6 
otherwise 



Since, by Lemma [3 l-O^I > 1/2 in -B(0,2) \ -B(0, 1), Dk > dk, and ^ follows from jMl) and the 
estimate > b. 

□ 
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Corollary 5. Assume that M > d. Then d! = d and dk = 1 for each k. 



Corollary 6. Assume that M < d. Then d' = M and d^ € 



d_ 

M 



d_ 

M 



+ 1 > for each k. 



Proof of Corollaries\^ and\^ By combining (|93p and (|94p we obtain the lower bound for Ff, in $7: 

M 



Fe{ve) + Ci > vr^ {deg5B(,^_5)(^;,)2| \n5\ + 62|degaB(a„5)(^e)|| InCl} • 



(100) 



i=l 



The conclusions of the above corollaries are obtained by solving the discrete minimization problem 
of the RHS of (HOO]). □ 

As a direct consequence of Proposition |3] and Lemma |8l we have 

Corollary 7. There is C > independent of e s.t. for 1 > p > 25 we have 



iVvel^dx = TT^dlllnpl +0(1) 



fc=l 

\d\ lnp\ + 0(1) in Case I 

M 

TT _ min d^l ln/)| + 0(1) in Case II 

_di,...,dM6Z 

di+,..+dM=d 



5.2 Existence of the limiting solution 

We now return to the proof of Theorems [1] and |2j 

Recall that {if , is a set of distinct elements of {1, ...,M}. We choose e small enough so 

that ij's are independent of e, thus we may simply denote this set by {ii, ...,id'}- In Case I, we have 
d' = d and we may assume that {ii, ...,id'} = {1, ■■■,d}. In Case II, we have d' = M. 

Lemma [7] and Corollary [7] imply that for an appropriate extraction e = En iO and for a compact 
K C Q \ {fflij, Oj^/}, there is Ck > s.t. for small e we have 



FeiVe,K)<CK 



and 



\ve{x)\ > I - C\ lne|~^/3 for all x e K. 

Therefore, when e —t- 0, up to a subsequence, there exists v* G H^(Q \ {oj^ , Oj^, }, S^) s.t. 
Ve^v* £ Hl^{n\{ai^,...,ai^,}). 

We now fix such sequence and a compact K G {flii, ...,aj^, }. U K G {a^, 1 < i < M}, 
then we have K Dujs = ^ ^or small e. By exactly the same argument as in Proposition [7] we deduce 
that Vs is bounded in C^{K) for all A; > and 1 — < Ck£^ in K. 

Consequently, up to subsequence we have for a compact set C \ {oi, om} 



V* in C^{K) and 1 - \ve\'^ < Cre^ 



(101) 



Now, assume that is s.t. K C il.\{ai^, fli^, } but K Dujs ^ (then we are in Case I). Without 
loss of generality, assume K = B{ako, -^)' where a^g G {orf+i, om} and i? > is sufficiently small 
in order to have K n {ai, om} = {ofco}- 
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Let /ig := tiQxVs- Since dK C \ {ai, ...,aM}, we have — t- /iq in C°°{dK) (possibly after 
passing to a subsequence). Since deg(/i£, = we have deg(/io, 9X) = and consequently there 
is some tpQ G C°°(aK,M) s.t. /iq = e*"^". 



Let i be a n.,™ of / |V„P the class //;„(/CS'). Cleaaly, 



K Jk 



On the other hand, since Ue < 1, we may construct (in the spirit of [TO] ) a test function and find 
that (see formula (93) in |10] ) 



F,{ve,K) < \ f \V^,\^ + Ce, (102) 



where ■i/'g is the solution of 



Here, ipf, is defined by 



As e — )• 0, we have 



AV'e = mK 
Ipe = on dK 



j^, = _^ on dK. 

\hJ 



Tpg — 7- V'o strongly in H^{K), where J ' . (103) 

[V'o = V^o on dK 

From the fact that in L'^{K), C/^ — t- 1 in Lp'{K) and \Us\ < 1 we have U^Vve v* in 
LP'{K). Consequently, we obtain 

-/ |Vv*p < liminf- / U'}\Vve\^ <\\mmi Fe{ve,K). (104) 
2 Jk 2 e->o 

Combining ffT02]) . ffT03]) and flT04]) we deduce that 



Ja' 



It follows that V* minimizes the Dirichlet functional in 

Hl^{K,S^) := {v G H^{K,S^),v = ho on dK}. 

We find that hence v = v* in K. By a classic result of Morrey [24] (see also [lOj), f* satisfies ([5]). 
Moreover, as follows from weak lower semicontinuity of Dirichlet integral, (|102p . (|103p and (|104p 

-/ I V?;* p < liminf- / \Vv,\^ < lira sup Fe{ve, K) <- [ \Vv*\^. 
2 Jk <^^o 2 Jj^ 2 

Therefore, 

converges to v* strongly in H^{K). (105) 

From (jlOip and (jlOSp we obtain that — t- v* in {{^^^{Q, \ {oi, ...,arf/}). The convergence up to 
dfl will be established in the next section. 

In order to prove Assertion 3. of Theorem [1] and Assertion 2. of Theorem [21 note that, for small 
p > 0, the estimate (jlOip implies that := trg5(a, p)^^ satisfies the conditions (|13p and (|14p of 
Theorem |3l This gives us 3. of Theorem [1] and ^. of Theorem [2] 

Assertion 5. of Theorem [2] is is a consequence of Corollary |6l 
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5.3 The macroscopic position of vortices minimizes the Bethuel-Brezis-Helein 
renormahzed energy 

Let us recall briefly the concept of the renormahzed energy Wg{{bi,di), {bk,dk)) with 

(geC^{dn,S^) s.t. degg^{g)=d 
61, bk £Q.,bi^ bj for j 
^di G 7, and Yli di = d 



For small p > 0, consider Clp = i}\ UiB{bi,p) and the minimization problem 

/p((6i, ...,6fc), (di, = inf ^/ |Vu;p 



w=g on dQ 



no 



Such problem is studied in detail in [llj (Chapter 1). In particular Bethuel, Brezis and Helein 
proved that for small /O, we have 

/p((6i, {di,...,dk)) = vrdi lnp| + di), {bk,dk)) + Op{l). 

This equality plays an important role in the study done in |1 1] . In the minimization problem of the 
classical Ginzburg-Landau functional 



2 

the vortices (with their degrees) of a minimizer tend to form (up to a subsequence) a minimal 
configuration for Wg. 

We prove in this section that the (macroscopic) location of the vorticity of minimizers of Fg, is 
related to the minimization problem of Wg((6i, 6^), (di, dk)) with 61, bk € {ai, oa/}- 
We present here the argument for Case I (Theorem [1]) . The argument in Case II is analogous. 
The proof of Assertion ^. relies on two lemmas, providing sharp upper and lower bounds. 

Lemma 9. There exists pQ > s.t., for every p < pQ and every e > 0, we have 

Feive) < TTd\ In p\ + dJ{e, p) + Wg{{a,,,l), (a^,, 1)) + Op{l), (106) 
where J{e,p) = inf„g^i^(5^(o)) Fe{u) with gp = ^ on dB{0,p). 

Proof. The proof, via construction of a test function, is the same as proof of Lemma VIII. 1 in 

o. □ 



Lemma 10. Let p > 0, p < pQ. Then for small e we have 

Feive) > nd\ lnp\ + dJ{e,p) + Wg{ia,,,l), (a^,, 1)) + Op(l). (107) 



Proof. Split the domain 0, into two sub-domains: \ UiB^ak^, p) and UiB^ak^, p)- We start with 
the lower bound in the first sub-domain. By the previous estimate, weakly converges to v* in 
H^{^ \ UiB{aki, p))- This implies that 



liminf- / U^\'^Ve\'^ > - [ \Vv*^^ 

2 JQ\Ut.B(a,. ,p) 2 Jq 



n\UkB(ai p) ^ Jn\UkBiai p) 



Here, we used the fact that, since f/^ — t- 1 in L (fl), {U^l < 1 and Vv^ Vv* in L (Q\[JkB{aif,, p)), 
we have UgVve Vv* in L'^{Q \ LlkB{ai^^, p)). 
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Thus we deduce that, for small e, 
1 



ln\UkB{aii^,p) ^ JQ\VJkB{ai^,,p) 

On the other hand, as proved m pT| . 

1 



Ue\^Ve\^>- |Vf*p-p2. (108) 



Vv*\' > ^dln - + Wg{ia,„l), (ai^l)) + Op(l). (109) 



/n\UfcB(a,^,p) P 

Thus, combining (jlOSp . (|109p and using Proposition EJ for e sufficiently small, we have 

1 



F,{ve, n \ Ufcfi(a,,, p)) > ^dln - + WgHai^l), (a^,, 1)) + Op(l). (110) 

By Theorem |4] and Corollary |5] we have the following energy expansion: 

F,ive,Biai,,p))=7rlnp + TTb''\lne\+TT{l-b'')\ln6\+W{a) + Wo{fo) + b^l + Oe{l). (Ill) 

Similarly, applying Theorem |4] to J{e,p) we obtain 

J{e,p) = TTln p + 7rb^\lne\ + 7r{l - b^)\ln6\ + W{a) + Wo{z/\z\) + b^-f + Oe{l). (112) 

Here, the local renormalized energy T^(q;) is given by (|90p and is the same in (jllip and (|112p . 

From ([73]), Wo(/o) > while Wo(]|) = 0. Consequently, we have Fe{ve, B{ai^, p)) - J{e,p) > 
Os{l)- Hence Vp > there exists Ep > s.t. for e < we have 

FeiVe,Bia^„p))>Ji£,p)-p^ 

and thus 

F,{ve,UkB{ai^,p)) > dJ{e,p)-dp\ (113) 

which gives the lower bound in the second sub-domain. From (jllOp and (|113p the bound (|107p 
follows. □ 

Combining Lemma |9] and Lemma [TOl we see that the points {aj^,l < k < d} minimize Wg 
among ai,...,aM- The expansion ([6| follows from (|106p . (|107p and (|112p . 

We next turn to convergence of up to the boundary. It suffices to prove the //^-convergence 
of Vs in i}p = Q \ L)mB{ai^, p) (for small p > 0). We argue by contradiction and we assume that 
there are some pi > and rj > s.t. 

liminf- / |Vu£p > - / \Vv*\'^ + ri. (114) 

Note that for all p < pi, (|114p still holds in ^Ip. 

If, in the proof of Lemma [TOl we replace (jlOSp by (|114p (with pi replaced by p), then we obtain 
for small p a contradiction with Lemma (9] The proof of Theorem [1] is complete. The last assertion 
of Theorem [2] is obtained along the same lines. 

A Proof of Proposition [2] 

Let xo G Vr be s.t. Br = B{xo, R) C 0,\lJs and assume that xq = 0. 
We follow the proof of Lemma 2 in jlO] . 
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In Bn, r] = 1 — Us satisfies 

-e^Ar/ + tr] = -'q{rf' - 3r/ + 2 - t) in 

r/ < 1 on SBr ' 

here, t will be cliosed later. 

Since r/ G (0, 1 — 6), if we take t = 6(1 + 6), then we have 

-e^Ar/ + t?? < in Br. 

On the other hand, the function w{x) = e^^^^^'^^^'^^ satisfies 

-e^Aw + tw= [-4e2-)/(l + 7|xp) + t] t(; in S/j 
w = 1 on 



A simple computation gives that 



Take 



-e^i^w + tw> {)m Br ^ ^ -e + Ve^TW 
7>0 ^0<^^ 2i^^i ■ 

-e + Ve^ + ti22 



> 0. 



Setting V = rj — w, we have 



-e^A?; + < in 5/? 



V < on 

By the maximum principle, we have f < in Br. Therefore, 

r/(0) < exp < — > <Ce ^ . 

Consequently, (|9]) holds in {x G | dist(x, dQ) > R, dist(x, cos) > R}. The estimate close to the 90 
is a direct consequence of < C/g < 1, (l9| holds in {x € O | dist(x,0r2) > i?, dist(2;, w^) > R} and 

the equation —AUs = —^Us(l — IC^P) in {x G $7 I distfXjCJ^) > R}. Using a similar argument, we 
establish (|9]) in the case VrDojs. The proof of <^ is complete. 

In order to prove (|10p . note that in Wr := {x G ^}\dist{x,dujs) > i?, dist(x, 957) > R} the 
function r] = Us — Us satisfies Aij = ^(a^ — [/|). Thus, applying Lemma A.l [10] to ry in conjunction 



with ([9]) and the fact that R> e, we obtain 



C ~~ 
|Vr/| < in Wr. 

Thus (|10p holds far away from dQ and the inclusions. 
We next prove that the bound (jlOp holds near dQ. 

Indeed, fix a smooth compact K C s.t. for small 6 we have uig C K. Clearly, by ([9|, 
< ??x := i^dKV ^ Ce~~. In \ IT, ry satisfies 

' Ar] = j^U{l + U)r] mn\K 
rj = on 00 

= on (9i^ 
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Let 7] = rji + 'i]2 he s.t. rn solves 



and r/2 satisfies 



Ar]i = ^Uil + U)r] mn\K 
r?i = on dn U dK 



Ar/2 = mn\K 
r/2 = on di} 
{V2 = i]K on dK 



Note that ||?/2||l°° < C'e e and thus < Ce e . 

Lemma A. 2 in |10) imphes the existence of a constant Cq\x > s.t. 



|V??i| < 



mn\K. 



In order to estimate Vr/2 near dfl, we express 7/2 in terms of Green's function G{x, y) in Q\ K: 
function, i.e. 

f dG 

mix) = - VK{y)^{x,y)dS{y). (115) 
JdK 



It follows from (|115p and (|9]) that |V7/2| ^ C^e e away from SiT. The estimate (jlOp is proved. 

B Proof of Proposition [3] 

This appendix is devoted to the proof of Proposition O 
We prove the first assertion: when M > d we have 

inf Fe{v,n) < vrdft^lln^l +7r(i|ln5| +0(1). 

Fix first d distinct points-centers of inclusions oi, a^. Let pQ := 10~^-min(dist(aj, 50), miuj^j \ai— 
aj\ > 0). Consider i; to be a smooth fixed function in \ Llf^iB{ai, po), such that \v\ = 1 in 
n \ utiB{ai,po) and 

{V = g on dil 

v{x) = ^\ ondB{a,,po)' 

Such a function clearly exists since the compatibility condition degQQ{g) = X^iLi degg^^^, ^^^^(t;) is 
satisfied. Let cq = 10^^ • dist(0, duj). For every 1 < i < M, consider a disc B{ai, cqS) C w^. By the 
choice of cq, we have dist((?a;5, i?(aj, cgf^)) > cqS. Therefore, using Proposition [2] 



-b^ < Ce-- in B{ai, cqS). 



(116) 



Consider the test function Vq defined as 



vl{x) 



V{X) 
X — tti 

\x — ai\ 

X — Qi 



for X € \ UiB{ai, po) 



for X e B{ai,po) \ B{ai,e) 
for x G B{ai,e) 



Using (|116p and ([H]) we have 

< rfjinel + 7rd(l-62)|ln5| + C = rf|ln^| +7r(i|ln5| + C. 



inf F,(t;,J]) < F,{vl) 
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Now we prove the second assertion: when M < d we have 

inf FJv,^}) < rflln^l +7r Vd^llnJI +C. 

Let di, dM e N be s.t. ^di = d. Set cq = IQ-^"^ • dist(0, du). For i G {1, M} s.t. di > 0, 
fix ai,i, ...,ad^^i £ B{0, W^cq) C w s.t. 



mm 



( min|a,-,i - a k, i\, dist {a j^i,du) I > Acq. 



Consider an e-dependent map Vq G H^{^ \ Ud^yoB{ai, 10°'co5),S^) s.t. 

Vq = g on dQ 

vf.(x) = ~ on dB(ai, W^coS) 

\x — Oip* 

and satisfying 



/C\Ud.>oS(ai,10'^co5) 

with C depending only on 0,a; and g. 

(Such maps do exist, e.g., consider the map introduced in Remark 1.5.) 



For i e {1, ...,M} s.t. di > 0, we consider a map vf € H'^{B{0, W^cq) \ ^%iB{aj^i,S,),S^) s.t. 
vl{x) = x'^VIxl* on dB{0, IO'^cq), 
vfix) = {x - aj^i)/\x - aj^il on dB{aj^i,S,), 

/ l^^ll^ ^ '^di \ In^l + C with C depending only on w. 



'B(0,10'*co)\u;iiB("j,>,0 

(For example, the map considered in Remark 1.5 in [llj has these properties). 

The necessary test function that satisfies the bound (|12p is obtained by rescaling the vf^s (in 
order to have maps defined in balls of size 5) and gluing the rescaled maps with Vq. 

C Proof of the ry-ellipticity Lemma 

The main argument in the proof of the ry-ellipticity result is the following convexity lemma which 
is a generalization of Lemma 8 in |9] . The proof of Lemma [11] is given in |15] . 

Lemma 11. [Convexity Lemma] 

Let C be a chord in the closed unit disc, C different from a diameter. Let S he the smallest of 
two regions enclosed by the chord and the boundary of the disc. 

Let O be a Lipschitz, bounded, connected domain and let g G C{dO,S). 

Assume that v minimizes Ginzburg-Landau type energy 

F{v) = j |d(a;)|Vi;|2 + /3(x)(l-|i;|2)2| dx 
in H^{0), with d,/3 G L°°(0,M) satisfying essinfd > 0,essinf/3 > 0. Then v{0) C S. 
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We prove the first part of tlie lemma [T] Let x G be s.t. dist(x, di}) > e^/^. We liave 
FM,B{x,s'/^))>^- [ \Nv,\' + ^il-\ve\'A 



'B(x,ei/4)\B{a;,ei/2) 
2 Jei/a r JdB{x,r) 

By Mean Value theorem, exists r E (e^''^,e^/^) s.t 



rf \\^v,\' + ^,{l-\vs\r}< 1,, , 

JdB{x,r) I £ J illnel 

There is C2 = C2(x,fe) > s.t if F^ive, B{x, e^^'^)) < lne|, we have 

Yai{vs,dB{x,r)) < C2X7 where Var (vg, 9i?(x, r)) := / l^rfel. (117) 

JdB{x,r) 

It follows that 

\ve\'^ >l-3C2Xon dB{x,r). (118) 

Indeed, arguing by contradiction, assume that there is e„ 4- ^-iid y„ G dB{x,r) s.t. |t'e„(yn)P < 
1 — 3C2X- Using (|117p we obtain that 

< 1 - C2X on dB{x,r) 

which implies that 



2(1-1) ^ 2_7rC7|r2;t2 
r 



2 

2\2 



< (l-l-.J^) 
JdB{x,r) 

^ lF,M.,B{x,en"')) ^ 8x\ 
l\lnen\ ~ b 

Clearly, the previous assertion gives contradiction. 

From (|117p and (jllSp . there is C = C{x,b) > and eo = eo(x) > s.t. for e < Eq, 

Ve:dBr^{zGBi\^z>l- Cx}. 

Using Convexity Lemma (Lemma 111 p . we find that \vs\ > 1 — Cx in B{x,r) D B{x,e^f'^). 
If dist(x, 89.) < e^/^, we denote Sr = 9r\dB{x, r), r e (e^/^, e-*^/^). Clearly, we have 

2 



b 



F,{v„B{x,e^/^)\B{x,ey^)) > 1^'^^^ ||V^;,|2 + 1(1 _ |,;,|2)2|. 



Using mean value argument and the facts that — )• 170 in C^(9r2,S^) and that < 1 — l^^l < e, 
there are r G {e^l^,e^/^) and Ci = Ci(||5o||ci> ^) s.t 

f ,2^ In I |2^2l ^ lF,{v,,B{x,e^l^)) + C^ 

r \\drVe\ +-J2{l-\V,\) — ■ . 

Using the same argument as before (taking = 90 B(x,r)) we obtain the desired result. 
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We prove the second part of the lemma. Let /Lt G (0,1) and x G {dist(a;, > £^/^}. Using 
mean value argument, there is r G {e^/"^ ,£^/^) s.t 



7 I 

JdB{x,r) I 



£ I 4 Ine 



4 

There exists Ci = Ci(/x,6) > s.t if £^4)) < Ci\ \ue\, we have 

\ax{ve-,dB{x,r)) < 1 — Iv^] < on dB{x,r). 

By Convexity Lemma l^gl > /x in B{x,r) D B{x,£^^'^). 

If dist(j;, < £V4^ denote Sr = nn dB{x, r), r G (eVs^ £1/4). since 

^F,(i;„B(x,£V4)\5(^^£i/2)) > r i.^dr/ /|Vt;,|2 + l(l-|t;,|2)2\ 

and using the conditions on g^, by mean value argument there is r G (£^/^,£^/^) s.t 

lF,{v,,B{x,ey^)) + C{\\go\\c^,^) 



r 



I [\drVe? + -^{l-\Ve\'f]< 

Jd(B(x,r)nn) I ^ ) 



ld{B{x,r)nn) I J 4Mn£| 

Using the same argument as before, the statement of the lemma follows. 
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